What the small angle CMB really tells us about the curvature of the Universe 
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It is well known that observations of the cosmic microwave background (CMB) are highly sensitive 
to the spatial curvature of the Universe, k. Here we find that what is in fact being tightly constrained 
by small angle fluctuations is spatial curvature near the surface of last scattering, and that if we allow 
A: to be a function of position, rather than taking a constant value everywhere, then considerable 
spatial curvature is permissible within our own locale. This result is of interest for the giant void 
models that attempt to explain the supernovae observations without Dark Energy. We find such 
voids to be compatible with the observed small angle CMB, but they must be either very deep (and 
unnaturally empty) or exist in a positively curved Universe. 



One of the great successes of modern cosmology 
has been the ability of Cosmic Microwave Background 
(CMB) anisotropies to constrain the geometry of the Uni- 
verse. A succession of ground, sub-orbital and space- 
based experiments [l| have led to increasingly tight con- 
straints on the curvature of space, k, when it is assumed 
to be a universal constant. However, in an inhomoge- 
neous universe k will not be constant everywhere, but 
will vary from place to place. Here we address the ques- 
tion of what CMB results imply if we allow k to vary with 
position. As a result, we place constraints on models of 
the Universe in which we live near the centre of a large 
under-density, or void. Observables in such models have 
been considered previously in , and have recently been 
used to explain the supernovae observations without re- 
course to Dark Energy [^-[l5l]. 

Of primary importance for constraining cosmological 
models are the CiS of the CMB angular power spectrum. 
These quantities are defined by an expansion in Legen- 
dre polynomials, Pe(x), of the form (<5T(n)<5T(n')) = 
h E<?( 2 ^ + l)C^(n-n'), where ST(n) is the CMB tem- 
perature anisotropy in the direction n, and angled brack- 
ets indicate an ensemble average. Here we will focus on 
the properties of the Cts on small angular scales. They 
are then a result of two processes: The imprint of cos- 
mological perturbations onto the last scattering surface, 
and the projection of that surface onto our sky. 

The first of these processes occurs early enough in the 
Universe's history that it is relatively insensitive to the 
effects of any spatial curvature. It can then be accu- 
rately described by linear perturbation theory about a 
flat background. The familiar set of peaks and troughs 
in the Cgs are then determined by cosmological parame- 
ters such as the expansion rate up to last scattering, and 
the relative densities of the different constituents of the 
Universe [l6| . 

The second process involves relating length scales at 
last scattering to angles on the sky today, and is highly 
sensitive to the geometry of the intervening space-time. 
Indeed, it is well known that non-zero k results in a shift 
of the acoustic power spectrum of small scale fluctuations 
in the CMB [17| , and that it is this effect that is responsi- 
ble for the stringent constraints on spatial curvature that 
usually imply k ~ 0. Such constraints, however, assume 



that k is a constant, throughout the Universe. Here we 
relax this condition, and allow k to vary with position, by 
considering the spherically symmetric Lcmaitre-Tolman- 
Bondi (LTB) space-time. We find that k(x) is only well 
constrained in the vicinity of the surface of last scatter- 
ing, and that even large local fluctuations in k will only 
produce moderate contributions to the shift. 

On small angular scales, the relative temperature of 
the CMB seen on an observer's sky in different directions, 
n, is given by 5T(£i) = A(nZ?Ls), where Dls is some 
measure of the distance to the last scattering surface and 
A is a solution of the Einstein-Boltzmann equations at 
the time of last scattering. In conformally static space- 
times, such as those with k =constant, Dls can unam- 
biguously be taken to be the conformal distance to last 
scattering r = sinh (y/—k J difj /y/—k, where dr\ = dt/a 
is conformal time. For more general space-times, how- 
ever, we will need to be more careful. 

Assuming that the radius of curvature is much greater 
than the scale of any perturbations, we have that the 
variance in temperature fluctuations is (ST(n)ST(h')) = 
J d 3 kVA(k)e ik < - & - & ') DLS , where we have defined the 
power spectrum to be -p A (|k|)(5 3 (k-k') = (A*(k)A(k')). 
What we arc ultimately interested in is the angle between 
vectors n and n' at the observer, d9 = cos _1 (n ■ n') ~ 
|n— n'|, hence, if dp = |n— h^D^s is the distance between 
two points at last scattering, then Dls — dp/ d6 = dA,LS, 
where dA,LS is the angular diameter distance to last scat- 
tering. This is a generic result valid for any curvature, 
constant or not. We will now approximate the C^s as 
a Fourier decomposition of (5T(n)ST(h')) over the sky. 
Defining the two dimensional wave number, q, such that 
q = |q| = i, we then have Cg ~ C q , where 

C q = [ dtf (ST(n)5T(n')) e^ 6 = -J—Va ( ■ 

J A,LS \dA,LsJ 

On small enough scales, of a few degrees and below, we 
expect this expression to be good enough for accurate 
parameter estimation [181 ]. 

Now consider two different space-times. Although it 
can be arranged that observers in each of these will wit- 
ness identical last scattering surfaces at identical red- 
shifts, the geometries between those observers and the 
surface will be different in each. Let us write dA,LS for 
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then gives the shift parameter as 



FIG. 1: Upper panel: The shift parameter, S, as a function of 
z\ for Q.k = 0.3 (solid), 0.5 (dotted) and 0.7 (dashed). Lower 
panel: S as a function of cosmic time, t\ = i(zi), for the same 
models (to is time today). 



the angular diameter distance in the first space-time, and 
Aa,ls for the angular diameter distance in the second. 
We can then relate the angular power spectrum in the 
first space-time, Ci, to that in the second space-time, 
Ce, via Ci = S 3 Ci/g, where S = (Ia.ls /d-A.LS is known 
as the shift parameter. This situation (of identical last 
scattering surfaces but different geometries) is often en- 
visaged when considering the effect of a non-zero, and 
constant, k. In that case both space-times are confor- 
mally static, and so it suffices to use the conformal (or 
optical) metric. The effect of k ^ is then to alter the 
conformal distance to the last scattering surface, and the 
Ct$> of each observer can be related by a shift parameter 
that is the ratio of these conformal distances . 

Now consider a toy model with a region of curved FRW 
extending out to some redshift, z\, in a universe that is 
otherwise flat. We have in Friedmann-Robertson- Walker 
(FRW) cosmology that the angular diameter distance is 
given by d,A = ar, where r is conformal distance (de- 
fined above) and a is the scale factor of the universe. A 
dust-filled FRW universe can be shown to have <Ia, as a 
function of the redshift z = ao /a — 1, given by 



d A = 



sinh(2f) 



(l + z)H y 



where f = sinh 



sinh 



(1) 



H is 



the Hubble rate, subscript denotes a quantity measured 
by the observer at z = 0, and f2fc = —k/a^H^. The shift 
in CMB peaks from a globally flat universe is now given 
by the ratio d A UIVcd /d A at at z\, when H has been matched 
at last scattering (and so is also matched at zi). In a flat 
universe we have H^^ &t = H 2 /(l + z) 3 , where H z is the 
value of H at redshift z, and in a spatially curved universe 
we have i? 2 curvod (l - Q k z/(1 + z)) = H 2 /(l + z) 3 . This 
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In Fig. [T]we plot this shift as a function of Z\ and the 
corresponding cosmic time for three choices of curvature, 
fifc = 0.3, 0.5 and 0.7. At large z\ we recover the famil- 
iar result that k < leads to S > 1, so that the acoustic 
peaks of the CMB are shifted to smaller angular scales. 
However, if we consider curved regions out to lower red- 
shift, then this result is no longer true: At Z\ % 4 nega- 
tive curvature causes S < 1. This is ultimately due to the 
presence of Ho in Eq. (Q]). Measuring d A in units of 
S would increase monotonically with z. The requirement 
that Hls is the same in both space-times, however, leads 
to different values of Ho in each. At low redshifts the ra- 
tio of these Hubble rates is great enough to cancel what 
would otherwise be a positive S — 1. 

In the lower panel of Fig. [1] we plot S all the way 
out to last scattering, this time as a function of cosmic 
time, t\ = t(zi), in the fiducial flat model. Here it can be 
seen that most of the shift parameter is due to geometri- 
cal effects shortly after the surface of last scattering, at 
z* ~ 1100, with any effects due to our local geometry 
contributing significantly less. In fact, for ~ 0.7 it 
can be seen that there is only a ~ 5% shift caused by all 
of the geometry out until the Universe was ~ 5% of its 
current age. The rest of the ~ 70% shift at last scattering 
is then primarily due to the geometry experienced by the 
CMB photons in the first ~ 5% of the Universe's history. 

Now let us consider models in which k is a smoothly 
varying function of position, as emerges in a universe with 
large density fluctuations 5] . To achieve this consider the 
LTB model, whose line-element is given by [20( 



ds 2 = -dt 2 



a|(t, r)dr 2 
1 — k(r)r 2 



a\{t,r)r 2 dn 2 



(3) 



where 02 = (rax)' , and primes denote partial derivatives 
with respect to r. The FRW scale factor, a, has now been 
replaced by two new scale factors, a\ and 0,2, describing 
expansion in the directions tangential and normal to sur- 
faces of spherical symmetry. These new scale factors are 
functions of cosmic time, t, and distance, r, from the cen- 
tre of symmetry, and obey a generalization of the usual 
Friedman equation such that 



8nG m(r) k{r) 



(4) 



where over-dots are partial derivatives with respect to t. 
The energy density is given here by p = (mr 3 )' /a,2a 2 r 2 , 
and redshifts by 1 + z = exp{/(d2/a2)<ii}, where the 
integral is along a past directed radial null geodesic. 

The LTB space-time is fully determined by a choice 
of the three free functions k(r), m(r) and to(r). The 
first two of these are specified above, and the third is the 
'bang-time', which in these models need not be the same 
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FIG. 2: Upper panel: SH (see text for definition) as a function 
of central void curvature, Q.k, for asymptotically flat voids 
with FWHM at redshift of 0.4 (solid), 0.6 (dotted) and 0.8 
(dashed). Lower panel: The shift parameter, S, as a function 
of central void curvature, Qk, for the same three voids. 

at all points in space. Without loss of generality, we can 
then make a gauge choice such that m =constant. We 
will also initially consider a simultaneous big bang, with 
to ^constant. These models have been much studied re- 
cently, as a space-time with local negative curvature al- 
lows for the possibility of explaining the supernovae data 
without Dark Energy. A fit to the data is often found to 
be a void with Qk ~ 0.7, and a width of z ~ 0.8. This is 
a significant amount of spatial curvature, extending out 
to large distances, and one may initially suspect that the 
sensitivity of the small angle CMB to spatial curvature 
may be sufficient to impose strong constraints on these 
models 

To investigate if this is indeed so, let us consider a 
negative local curvature perturbation in an otherwise flat 
universe. An observer at the centre of such a void will 
see a last scattering surface at z*, and can straightfor- 
wardly calculate H at this surface in terms of their lo- 
cally measured value. We also require a fiducial observer 
in an FRW universe who will witness an identical last 
scattering surface, with the same Hls- To ensure that 
these observers use comparable measures of distance we 
will enforce the conditions that they have the same lo- 
cal geometry [2l|. This choice ensures that distances to 
nearby co-moving objects are the same. We also require 
that they both see last scattering surfaces at the same 
so that any effects due to the redshifting of solid angle, 
for example, are automatically included. 

The shift between the open FRW universe and the void 
model is then given by S\ — ^a T ls/^TlS' wnere the 
angular diameter distance in LTB is given by d^^ s = 
a-i,LS r LS, and in the FRW universe by Eq. ([I]). In the 
case of the void model, the values of a l5 r and H at last 
scattering are found by integrating a radial null geodesic 
out to Zjf, using the solutions to Eq. (UJ). Hq in the 
open FRW universe is then found by taking the same 
Hubble rate at last scattering as in the LTB model, and 
propagating it forward until today in the FRW geometry. 
Of course, we know the shift parameter between open and 
flat FRW universes, S 2 = d°l™ s / from ©, and so 
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FIG. 3: The 68% and 95% confidence regions from I > 100 on 
the central curvature, Qk, and FWHM of k in redshift, ze, of 
an asymptotically flat void with simultaneous big bang. Shell 
crossing occurs in the excluded hatched region. 



we can calculate the acoustic spectrum witnessed by the 
observer in the void in terms of a shift, S = S1S2 = 
^als/^a^lS' f rom a spatially flat FRW model, and a 
change in Hubble rate, SH = Hg ds /Hfr TB . 

The shift, S, and change in Hubble rate, SH, for an 
asymptotically flat void formed from a negative Gaussian 
perturbation in k(r), are shown in Fig. O We find that 
a good fit to the WMAP data exists for S ~ 0.9 and 
SH ^0.5, and so a void model will need to be capable of 
achieving similar values if it is to be considered viable. It 
can immediately be seen that for moderately deep voids, 
with Qk ^5 0.9 at the centre, both S and SH deviate 
insufficiently from 1 (22J. It can also be seen that S is 
largely insensitive to the width of the void. In light of 
what we considered above, these results can be easily 
understood: Most of the contribution to the shift does 
not come from the local geometry, but from early times 
when the CMB photons were well outside the void. The 
discussion above also explains why the presence of the 
void shifts the acoustic peaks to larger scales, rather than 
smaller. 

We may, however, wish to consider more extreme voids 
in which we allow Q). > 0.9 at the centre. In this case 
both S and SH can deviate considerably from 1. We have 
used a Monte Carlo Markov Chain algorithm running a 
modified version of the publicly available CosmoMC [23| 
to find constraints on Qk an d the extent of the void, Ze 
from the WMAP 5 year data, with I > 100. Ignoring 
models with shell crossing [24j we find the best fit void 
has Qk = 0.986+° °oi a t its centre, when we consider voids 
with a FWHM in k(r) at less than z = 1.2. The likeli- 
hood contours in this parameter space are shown in Fig. 
[3J The best fit model in this region of parameter space is 
effectively identical to ACDM, with Ax 2 = 0.4. This suc- 
cess at fitting the detailed shape of the CMB adds some 
perspective to claims that the CMB tightly constrains the 
geometry of the Universe and provides evidence for Dark 
Energy (when accompanied by an independent measure 
of Hq). In performing this fit we have used the Hubble 
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FIG. 4: Top panel: the Ces for the void in a non-flat back- 
ground (dashed) and for ACDM (solid) are essentially indis- 
tinguishable. Bottom panel: the distance modulus for the 
same two models. Data points are from the WMAP 5 year 
data [lj (top panel) and SNLS [27| (bottom panel); in the 
latter case, the position of the data points move as they are 
fitted to the two models (see [3). 

key project constraint h — 0.72 ± 0.08 [251 ]. The fit to 
the large I CMB then gives us the combined results for 
local Hubble rate h = 0.70 ± 0.01, and baryon fraction 
il^h 2 = 0.07+Q Such extreme voids, however, will only 
fit the supernovae data if their shape is fine tuned. 

Of course, one will be interested in more general void 
models. In particular, it is possible to conceive of a 
void in a spatially curved FRW universe, instead of a 
flat one. In this case one is subject to the familiar sen- 
sitivity of the CMB to spatial curvature, and we have 
verified that S can effectively be set to any value with a 
suitable choice of asymptotic curvature [26(. The value 
of SH, however, is not so sensitive to k in the background 
space-time. To achieve SH ~ 0.5, with a simultaneous 
big bang, one still needs, then, an extremely deep void. 
Alternatively, if one is prepared to abandon the notion 
of a big bang happening at the same time at all points 
in space, then one can straightforwardly achieve a larger 
contrast between local and asymptotic Hubble rates. To 
this end, we find that a Gaussian void embedded in a 
closed universe of Slfc = —0.05, that has a FWHM in k 
at z — 0.31, in to at z — 0.92, and with an age of the 



universe in the centre of the void that is 87% of that in 
the asymptotic regions, we can fit the Cis just as well 
as the best fit model above, with a much more modest 
fij? cal = 0.61. The CMB acoustic spectrum and distance 
modulus plot for this void are shown in Fig. [4] together 
with the ACDM best fits. Changing the detailed shape 
of the under-density will change the numbers involved 
above, and, in particular, if one can find other voids that 
allow SH ~ 0.5 then these models will very likely provide 
a good fit to the data too (with the appropriate choice of 
background curvature, to give the correct shift). 

In conclusion, we find that the observed acoustic spec- 
trum of small angle CMB fluctuations is primarily only 
sensitive to curvature at high redshifts. Local curvature 
has much smaller, and even opposite, effects. By consid- 
ering LTB models, in which k — fc(r), we demonstrate 
that large local fluctuations in spatial curvature produce 
only moderate shifts in the CMB acoustic spectrum. As 
a result, the local void models that seek to explain cos- 
mological observations without Dark Energy are not au- 
tomatically ruled out. Fitting to the WMAP 5 year data 
shows, however, that the simplest voids (with simulta- 
neous big bang, and asymptotically flat) are required to 
be extremely deep, with flk ~ 0.99 at their centre. Such 
depth is needed to achieve an admissible Hq, and to pro- 
duce a sufficient shift in the location of the spectral peaks. 
Alternatively, we can give up asymptotic flatness and a 
simultaneous big bang. By embedding the void in a suit- 
ably curved background it is then possible to shift the 
acoustic spectrum by any amount. Similarly, the local 
Hubble rate is sensitive to the local bang time, and by 
making our region of the universe younger we can in- 
crease Ho. There then exist some void models that fit 
the WMAP 5 year data just as well as ACDM. However, 
if we really do live in a large, local under-density in the 
Universe, it may have to be more complex than previ- 
ously thought in order to account for CMB observations, 
as well as supernovae. 
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